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Abstract 

We study a simple model of a quantum Hall system with the elec- 
trons confined to a linear, narrow channel. The system is mapped to 
a ID system which in the low-energy approximation has the form of a 
Luttinger liquid with different interactions between particles of equal 
and of opposite chiralities. In a previous paper (Part I) we studied this 
mapping at the microscopic level, and discussed the relation between 
the parameters of the 2D system and the corresponding ID Luttinger 
liquid parameters. We follow up this study here and derive the low- 
energy form of the 2D electron correlation function by applying the 
mapping from the ID correlation function. We examine in particular 
the modification of the asymptotic behaviour of the correlation func- 
tion due to interactions and find for the 2D function a similar fall off 
with distance as in ID. 

1 Introduction 

In Part I [1] we have examined in detail the relation between the quantum 
Hall system and the Luttinger liquid in a simplified case, where an explicit 
microscopic mapping between the two systems can be performed. This is a 
case of integer filling where the 2D electron gas is constrained to a narrow 
channel by a harmonic oscillator potential, and where the electron interac- 
tion is assumed to be sufficiently long range to interconnect the two edges. 
The corresponding ID system is in the low energy approximation described 
by a full Luttinger liquid with both chiralities present, and with different 
interactions between particles of the same chirality and with opposite chiral- 
ities. In [1] we in particular focussed on how the Luttinger liquid parameters 
were modified by the electron interaction. 
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In the present paper we study the electron correlation function and show 
that an explicit form of the 2D correlation function can be found in the low 
energy approximation. As is well known the electron correlation function 
of the ID Luttinger liquid can be derived by the bosonization technique 
[2, 3, 4] and we follow [4] in order to establish this for the present case. The 
mapping discussed in [1] is then used to derive the 2D correlation function. 
The function is found to have a natural separation into two parts, which we 
associate with the bulk and the edge of the 2D electron gas. 

We focus in particular on the effects of the electron interaction on the 
correlation function. A part of the motivation for this is the recent discus- 
sion about the effects of interactions on the electron correlation function in 
the fractional quantum Hall effects, where numerical studies [5, 6, 7] have 
indicated, in some cases, deviations in the asymptotic behaviour of the cal- 
culated 2D correlation function from the expected behaviour based on the 
correspondence with the ID chiral Luttinger liquid [8]. These studies in 
turn have been motivated by the difference between the observed tunneling 
resistance and the expected resistance based on the chiral Luttinger descrip- 
tion of the quantum Hall edge [9, 10]. The discussion illustrates the point 
that the analysis of the edge effects to some extent has to rely on qualita- 
tive arguments and on numerical studies of systems with a limited number 
of particles, and that this may leave room for some uncertainties due to 
different interpretations. 

The simplified model studied in this paper cannot be directly related 
to the unsettled questions concerning the edge effects at fractional filling. 
Even so, it seems interesting to examine in detail a case where an explicit 
mapping between the 2D quantum Hall system and a ID Luttinger liquid 
can be performed, and to examine in this case how the electron correlation 
function of the quantum Hall system is related to the correlation function of 
the Luttinger liquid. As a particular result we find no difference between the 
asymptotic parameters in ID and 2D, caused by the interactions. However, 
the region where deviation from the asymptotic form is non-negligible, may 
extend to a distance that is rather large when measured in units of magnetic 
length. 

As discussed in [1] we make certain simplifying assumptions for the model 
studied. We assume a soft harmonic confining potential, with an oscillator 
frequency that is small compared to the cyclotron frequency, and also as- 
sume a weak particle interaction in the form of a gaussian potential. In the 
following we first make a brief review of this model and its mapping from 
2D to ID. We further discuss the electron correlation in the ID description 
and show how the 2D correlation function can be derived. We supplement 
these results with some graphs that demonstrate how the electron correla- 
tion function and the density function in 2D is modified by the interactions 
for some specific choices of parameter values. 
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2 The model 



The single particle part of the 2D electron Hamiltonian has the form 

H = ^(P* + eBy) 2 + ^Pl + \™y (1) 

which describes spinless electrons in a constant perpendicular magnetic field 
B and subject to a harmonic confining potential in the y-direction. Landau 
gauge has been chosen here, with vector potential 

A x = -yB , A y = (2) 

We re- write this in terms of the effective cyclotron frequency u) c = \/ lo 2 + w2 > 
with lo c = eB/m, 

H = 2^ p y + 2 m ^ {y + ^! Px) + 2^ (3) 

The condition for a soft confining potential, uj « uo c is assumed in the 
following. 

Since p x and H commute they can simultaneously be quantized, and 
also the operator defined by the first two terms of (3) can be quantized 
simultaneous with p x . This part of H has the form of the Hamiltonian 
for an electron in the effective magnetic field B = \J B 2 + m 2 u) 2 /e 2 . As 
a basic assumption we assume throughout the paper that the low energy 
approximation is well satisfied, where electrons are confined to the lowest 
Landau level of the Hamiltonian for the effective magnetic field. This means 
that a complete set of eigenstates can be written in the form, 

Mx, y) = -^=e ikx Mv - Vk) (4) 

with hk as the quantized value of p x , L as a normalization length in the 
x-direction and yu = —{uj c /uj c yP B k as a k-dependent shift in the y-direction. 
V>o is the ground state wave function for the (effective) harmonic oscillator 
equation in the y-direction, 

My) = (^e'^ y2 (5) 

In the expressions above the magnetic length of the effective magnetic field 
has been introduced, Ib = \J fr/eB 

The restriction to the lowest Landau level effectively reduces the dimen- 
sion of the system from 2 to 1, which is demonstrated by the fact that 
the wave functions (4) are characterized by a single momentum variable k. 
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An explicit mapping from two-dimensional wave functions ip(x,y) to one- 
dimensional wave functions can be written as 



ip(£) = / dxdyf(£-x,y)ip(x,y) 




(6) 



with the integration kernel defined as 



f{x -£,y) = (x, y\0 = jJ2 e ik(X -®My ~ Vk) 



(7) 



k 



In the limit L — > oo the exact expression for this function is 



f{x,y) 



A 




(8) 



exp 



with 



A 




(9) 



We note that A is effectively a scaling factor in the x-direction, and as shown 
by the expression (3) such a scaling is indeed needed in order for the two 
first terms of the Hamiltonian to have precisely the form of the Hamiltonian 
for a charged particle in the effective magnetic field B. 

The mapping between wave functions can be used to derive the corre- 
sponding mapping between operators in the 2D and ID description. We 
give the expression only for the mapping of a local density interaction in 2D 
and refer to [1] for more details. Thus, the interaction in 2D we assume to 
have antisymmetric matrix elements in the coordinate representation of the 
form, 



F(n,r 2 ;riy 2 ) = V{v l - r 2 )(«5(n - r[)S(r 2 - r' 2 ) - 5(n - r' 2 )<5(r 2 - r' 1 ))(10) 



which means that it is defined by a single function of two variables, V(x, y) = 
V(r\ — r 2 ) with (x,y) as the cartesian components of the relative position 
vector ri — r 2 of two electrons in the plane. This operator is mapped into a 
non-local operator in ID, of the form 



where £ and £' represent relative coordinates of the two particles in ID. 
The function <?(£,£'; is determined by the one-particle transformation 




(11) 
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function (8), and has the form 
A 2 

g(tt';x,y) = - Aexp 
(V2vr l B ) 6 



._L {A 2 (x _ ^2 + A 2 (x _ + 2?A(e _ ^ )y)} 



with A denoting antisymmetrization with respect to the coordinate trans- 
formations £ — > — £ and £' — > — 

As a special case, used in [1], and applied also here, we consider the 
gaussian interaction 

V(v) = V e- a2r2 (13) 

The corresponding ID interaction has also a gaussian form, but depends 
on two rather than one length scale. In the momentum representation the 
interaction has the form 

V(k, k') = V ^exp[-^(fe - k'f - b 2 (k + k') 2 \ (14) 

where k and k! are the variables in the Fourier transform corresponding to 
£ and while Vq, a and b are constants given by 



V = 2 
a = A 



7T Vq 



l + 2a 2 P B a 
a 



^A 2 + 2a% 



b = °% (15) 



3 Low energy Hamiltonian and correlation func- 
tion in ID 

In the ID representation the second quantized Hamiltonian has the form 

h 2 1 1 
H = 12(^M k2 + 2 toc ) c fe Cfc + 4L S K ( fc ' fc/ )4i c fc 2 c *a-9Cfc 1 + g (16) 

A; q,ki,k2 

with k = (ki — k<i)j2 and fc' = (&i — ^)/2 + with M = m uj 2 /uj 2 as the 
effective ID mass and with V(fc, fc') given by (14). The operators ct and c& 
denote in the usual way electron creation and annihilation operators. As- 
suming a weakly interacting system, a low-energy form of this Hamiltonian 
can be found. Thus, with the assumption that excitations are energetically 
restricted to states close to the Fermi points ±kp, i.e., to states with k- 
values restricted by \k ± fejr| << kp, the non-local interaction (14) can be 



(12) 
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represented by two different local density interactions, V±(q) between parti- 
cles at the same Fermi point and V^g) between particles at opposite Fermi 
points. For the gaussian interaction these are given by [1] 

V 1 (q) = F exp[-£^] , V 2 (q) = exp[-46 2 4] V 1 (q) (17) 

If furthermore the first term of (16) is linearized in k around the Fermi 
points, the Hamiltonian can be brought into the Luttinger form 



H = v F h^( X k ~ k F ) : c\ k c x , k : +— ^{Vi{q)p x , q p x ,- q + V 2 {q)p x , q P- x ,- q ) 

(18) 



4L 



where the c^c operator and the particle density operators p x ^ q are normal 
ordered relative to the non-interacting ground state. In this expression \ 
denotes the chiral quantum number which initially is defined with % = +1 
for positive k and \ = — 1 for negative k. However, this relation between 
k and x can subsequently be relaxed, so that in (18) k runs from — oo to 
+cx) for both chiralities, without the low energy part of the theory being 
affected. Furthermore, in the linearized kinetic term the Fermi velocity has 
been introduced in the form 

vf = v f - ^-Al-e-^) (19) 

where v F = k F /M is the unperturbed Fermi velocity and where the interaction- 
dependent contribution to v F arises from the modification of the background 
potential due to presence of the Fermi sea [1]. 

The bosonization technique can be used to diagonalize the Hamiltonian 
(18) and to find the electron correlation function of the interacting ground 
state. As stressed in [4] one needs to consider the q = component of the 
charge density separately and to identify the q ^ of the charge density 
with the bosonic operators. In our notation the relations are 




where a q and a q satisfy standard boson commutation relations. 

Expressed in terms of the bosonic operators the Hamiltonian takes the 
form 

H = ^(v N (N-N ) 2 + vjJ 2 ) 



+ 



h 

2 



^ ^ VF q<lq aq q > ~4^h^ q q aqa ~ q > 



(21) 



6 



where unessential terms (constants or terms proportional to the particle 
number) have been subtracted. In the above expression we have included 
the particle number iV and the chiral charge J, 

N = N + J2px,o, J = Y.*Px^ ( 22 ) 

X X 

with Nq as the particle number of the ground state. The two new velocity 
parameters that are introduced are 

v N = v F + p-(l + e-^) = v F + p- h e-^ 

vj = v F + ^(l-e- Ab2k2 F) = v F (23) 

where only is modified by the interactions, as discussed in [1]. 

The Hamiltonian (21) is in the standard way diagonalized by a Bogoli- 
ubov transformation, which transforms the annihilation operators as 

a q = S^bgS = cosh Q q b q + sinh ( q b ] _ q (24) 

with 

S = exp[^ ^2 ( q {a\a)_ q - a q a^ q )} (25) 

q^O 

and C q defined by 

tanh(2C 9 ) = - - ^ - = -e-^l % (26) 
Expressed in terms of the transformed bosonic operators the Hamiltonian is 

q^O 

where again a constant has been subtracted. The Hamiltonian has a free 
field form, with the boson frequency given by 

u q = V(Vi(q)+^frv F )2-V 2 (q)2 \q\ (28) 

The bosonization technique makes it possible to determine the electron 
correlation of the ID interacting system. By use of the mapping between 
the ID and 2D system it also determines the correlation functions of the 
quantum Hall state. We examine this in the following and in particular 
compare the asymptotic behaviour of the correlation function in ID and 
2D. We also examine the effect of the interaction on the particle density. 
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The chiral (equal-time) electron correlation function of the ID system is 

c x (e-eo = (oi^t ( ^ x (05t|o) (29) 

where ip x ( x ) 1S ^ ne chiral electron annihilation operator, |0) is the ground 
state of the non-interacting Hamiltonian and S gives the transformation be- 
tween ground states of the non-interacting and the interacting Hamiltonian. 

Based on the bosonic form of the electron operators and the commutation 
relations between the bosonic creation and annihilation operators and S, an 
explicit expression for the correlation function in the interacting case can 
be found. We refer to [4] for more details (note however a difference in sign 
convention), and only cite the expression for the chiral correlation function, 

C x (£ - O = - ——, e" ) (30) 

with e — > + and where the function A{£) is determined by the interaction 
potential as 



Mi) = 4^ r sinh 2 C gS in 2 (^/2) (31) 



The asymptotic behaviour of the correlation function, for large £, can be 
found by considering the derivative of the function A(g), which in the limit 

L — > oo is 



dA(Q 



poo 

2 dq sinh 2 ( q sin(^) (32) 
■/ 



For large £ the integral collects the main contribution from small q, and we 
get the asymptotic approximation 

with 

Vb + 4vr^F , , 

7 = , (34) 

^(Vb + 4vr^ F ) 2 - exp(-86 2 4)y 2 
The corresponding asymptotic form of the chiral correlation function is 

C x (0 « r 7 (35) 

with K as a constant. The expression shows that the deviation from the 
non-interacting value 7 = 1 is determined by the strength of the interaction 
between the two Fermi points, which means between the two edges of the 
quantum Hall channel in 2D. 
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The full electron operator can be expressed in terms of the chiral oper- 
ators when the fourier components for non-negative k values are identified 
with the corresponding positive chirality components and the fourier com- 
ponents for negative k are identified with negative chirality components 

c k = c +j k if A; > 

c k = c_ fc if k < (36) 

This gives the following transformation between the chiral and non-chiral 
electron wave functions 

^)=E [ 2 L F x (C)^ + C)dC (37) 

with 

X e i L < > x 1 

Fx(C) = 2iL sin(f C) " ixe ~ ^ C^ix~e (38) 

where the last expression corresponds to the L — > oo limit. Since the F x 
function satisfies the following property 

F X (Q = J d£F x (Z + ±C)^(£ - l -Q (39) 

the corresponding mapping between the chiral and non-chiral electron cor- 
relation functions can be expressed as 



C M0 = E/ d(F*(()c x (t; + 



x 

Elf 1 p -*xM£+C) , 

A / t^— * . , e- A ^ +0 (40) 

x 4vr 2 7 C + «Xe e + C + «Xe' 

To evaluate this integral we express the factor e~ A ^ + ^ in terms of its Fourier 
components and perform the £ integral for each Fourier component by clos- 
ing the integration contour either in the upper or the lower half plane, de- 
pending on the sign of the ^-dependent exponent in the integrand. The 
result can be written as 

wo-En^-TT-)'-* (41) 

where the factor e~ A ^ is derived from e~ A ^ by canceling the contributions 
from the high frequency Fourier components for k < —kF (or k > fop). 
Assuming these high frequency components to be unessential we have A(£) ~ 
A(£) and the ID electron correlation function gets the form 

cmo = E c x(o - m = i^M e - A «) (42) 

x 

in accordance with the corresponding expression in [4]. 
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4 The electron correlation function in 2D 



We focus now on the problem of finding the electron correlation function 
for the 2D quantum Hall system. This function can be determined from the 
ID electron correlation function by use of the explicit mapping that exists 
between the ID to 2D systems. 

Before proceeding we introduce two simplifications. The first one is to 
set A = 1. This can either can be seen as an approximation which is valid 
when uj « uj c , or it can be viewed as a change to rescaled variables in 
the x-direction, where A is absorbed in the the new coordinates x and £. 
(However, we shall use the same notations for these variables as before). 
Since in the following we only refer to the magnetic length of the effective 
magnetic field, we also introduce the simplification of referring to this length 
as Ib rather than Ib- 

The correlation function of the quantum Hall state, which is translation- 
ally invariant in the x-direction, we write as 

C 2 d(x; y, rf) = U\x, y + r//2)^(0, y - rj/2)\ (43) 



This is related to the ID correlation function by a transformation of the 
form 

C 2 d(x; y,v) = J <K£; V, v) C 1D (x + £) (44) 
where the transition function <f> is determined from the /-functions as 
<H£; V,V) = J dx f*(x - £/2, y + rj/2) f{x + £/2, y - r,/2) 

' exp[--i-(^ + 4^)+r? 2 )] (45) 



" 2vr/| HL All 
This gives the following integral expression for the correlation function 

C 2 n(x;y, V ) = 2^/* exp[-^(^ + 4iy) +^)]™^^ c -^(46) 

With the interaction dependent function A(£) given by (31) this inte- 
gral cannot be evaluated explicitly. However, we note the presence of the 
gaussian function that effectively restricts the integration to an interval of 
length Ib- If the integrand is rewritten as the product of this gaussian func- 
tion and a function that varies slowly over the length scale Ib, the latter can 
be expanded about the centre of the gaussian to give a rapidly converging 
series where the lowest order terms can be evaluated. Since the function 
sin[fcp(x + £)] is rapidly oscillating, it is useful to rewrite the integrand in 
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terms of the two chiral contributions, 

C 2D {x;y,v) = ~ 4 l U 2 ^Zx&M-^tV 2 - ixkFx] 

B x B 

xjdt exp[-^-(£ 2 + 4i& + 4i X M|fl] g + g ± - £ } (47) 

where a shift ±ie (e — > + ) has been introduced for the location of the poles 
relative to the real axis. Note that the sign of this shift can be chosen freely, 
but has to be the same for both chiralities in the sum. 

We introduce a new, x _ dependent, variable in the y-direction, y x = 
V + xk F l%- F° r X = +1 it measures y relative to the edge located at —k F l 2 B 
and for x = — 1 it measures y relative to the other edge at +k F l 2 B . We also 
introduce a new, x-dependent, integration variable £' = ^ + 2iy x . This gives 
the following expression for the correlation function, 

C 2D (x;y,rj) = ~ 4 2 ^2 Y x{exp[-^-(?? 2 + 4j/%) - ix**^] 

B B 



Sc. 



1 „ p ->l(x+£'-2ij/ x ) 

< exp[ t^g 2 }——, } (48) 



with C x indicating an integration path shifted by 2iy x relative to the real 
axis. This integration path can be shifted back to the real axis, but due 
to the crossing of poles through the integration path during the shift, there 
will be compensating pole contributions. We therefore write the result as 

C 2 d(x; y, rf) = C int (x; y, rj) + C po i e (x; y, rf) (49) 

with the integral contribution given by integrals along the real axis, 

C in t(x; y, r]) = - 47t2U 2 Yl x{exp[-^2-(r? 2 + 4y 2 ) - i X k F x] 

-A(x+£-2iy x ) 



r 1 „ p-A(x+£-2iy x ) 

j^^-W^l^-} (50) 



while the pole contribution is found to be 

C po ie{x;y,ri) = 



y > k F l\ 

\r exp[- i ^ 2 -(x 2 + rf 2 - Aixy)} - k F l% <y< k F l% (51) 



y < -k F l 2 



B 



We note that there is no pole contribution for y > k F l 2 B and for y < 
—k F l 2 B . This is most easily demonstrated in (48) by choosing the shift 
— ie in (48) for y > k F l B and +ie for y < —k F l B . The integration paths 
then can be shifted to the real axis without crossing the poles. In the case 
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Figure 1: The integration paths for three different cases, a) y lies between the 
edge values ±kplg, b) y is smaller than — kpl\ and c) y is greater than +kpl%. C + 
denotes the integration path for the positive chirality term and C_ for the negative 
chirality term, and the location of the poles of the integrands is indicated for each 
of these by the nearby star. In case a) the integration paths, when shifted to the 
real axis (vertical arrows) , have to be compensated for by a pole contribution from 
the pole close to C+, while in the cases b) and c) the paths can be shifted to the 
real axis without passing the the associated poles. 



—kpl 2 B < y < kpl B a crossing of one of the poles cannot be avoided, and 
that gives rise to the pole contribution (see Fig. 1). 

The pole contribution C po i e , which is present only for values of y within 
the two edges, can in fact be interpreted as the bulk contribution to the 
correlation function, since it is independent of the location of the edges. 
The edge contributions are then given by Ci n p As shown by (51) the bulk 
contribution is insensitive to the electron interaction. 

The discrete change in C po i e at y = ±kpl B may suggest that the full cor- 
relation function C<id also changes discontinuously at the edges. However, 
that is not a correct conclusion, since the discontinuity of C po i e only com- 
pensates for the discontinuous change in C{ n t when one of the poles passes 
the integration path when y is changed. Thus, C2D(x;y,r/) is everywhere a 
continuous function of y. 

We consider next the contribution Ci n t and note that it essentially con- 
sists of two independent contributions, corresponding to the two chiralities 
X = dbl. Thus, each of these contributions is exponentially damped when y 
moves away from the corresponding edges at y = ^kpl B , and with Ib as the 
damping length there is little overlap between the two contributions when 
the separation between the two edges, W = 2kpl B , is much larger than Ib- 
In a similar way the integration variable is effectively limited to £ < 21b by 
the gaussian factor, so that both y x and £ are restricted to be of order Ib- 

The function A(x), which depends on the particle interaction, is a slowly 
varying function of x, with the interaction length 1/a as the typical scale on 
which the function changes. Since this is assumed to be much longer than 
Ib, an expansion of A(x + £ — 2iy x ) in powers of £ — 2iy x should give rise to 
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a rapidly converging series. Here we restrict the expansion to second order, 

A(x + £ - 2iy x ) = A(x) + a(£ - 2iy x ) J B(x) + ^a 2 (£ - 2iy x ) 2 C{x) + ... (52) 

where the higher derivatives of A(x) are expressed through the dimensionless 
functions B(x),C(x), .... This is essentially an expansion in powers of the 
small quantity al B . 

When the expansion is restricted to second order, the exponent of the 
integrand in (50) is quadratic in the the integration variable £, like in the 
non-interacting case, and the the effect of the interaction can then be ab- 
sorbed in the x and y dependent coefficients, 

C int {x;y,rf) = ~^rjT ^2 X {exp[-^2"(r/ 2 + 4y 2 ) _ i x k F x - A{x)\ 

B x B 

eM~M 2 ~ " 2i Vx)B(x) - \a 2 (i - 2iy x ) 2 C{x)} 

d£ =2 (53) 

x + i-2iy x V ' 



I 



By introducing a shift in the variable £ the integral can be written in the 
form 



cw*; ».i) = -4^E^* dt —z (54) 



Ex<' [ 'l <"f " 

*B X 

with coefficients 

$ = --^ T {r ] 2 + 4yl)-A(x) + 2a 2 y x C(x) + a 2 l 2 B B(x) 2 + ixkFX + 2iay x B(x) 

z = -^{l + a 2 l 2 B C{x))-i^{l-a 2 l 2 B C{x))-l B aB{x) (55) 
21b Ib 

where terms of higher than second order in al B are neglected. 

The integral in (54) can be expressed in terms of the error function in 
the following way, 



/ 



e _ 2 

dt ^— j- — = — iire z (erf (iz) + sgn(Im(^)) (56) 



where the term proportional to sgn(Im(z)) has a discontinuous y-dependence 
that exactly matches the omitted term C po i e . To demonstrate this, we note 
that the corresponding contribution to Cm can be written as 



z2 sgn(Im(z)) = -— ^ Yl Xexp[-^-(x 2 + rj 2 - 4ixy)]sgn(y x ) 

B x B 

C po ie(x;y,ri) (57) 



^ X B X 
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To obtain the correct exponent on the right hand side of this equation one 
has to evaluate <!> — z 2 to second order in oIb, consistent with the approx- 
imation we apply, and to note that x should be treated as of order l B , due 
to the exponential damping factor. 

For the contribution to Cj nt that is proportional to erf (iz) the situation 
is different. There is no exponential damping of this term with increasing x 
since the factor exp(— x 2 /l 2 B ) is compensated by a corresponding exponential 
increase of the error function. The overall dependence of this term with x 
is therefore a non-exponential decrease. 

For the full electron correlation function we now get the following ex- 
pression, 

C 2 D(x;y,r ] ) = -^{etfii z + (x,y)\ - erf[i z.(x, y)]} e - A ^) (58) 

with 

z±(x,y) = (1 + a%C(x))^- - l B aB(x) - i(l - a 2 l 2 B C{x)) V - M * (59) 

Mb Ib 

and 

2 2 

A(x; y, 77) = (1 + 2a 2 l 2 B C(x))^ T + A(x) - aB{x)x + -JL. - ^xy (60) 

This final expression for the correlation function includes effects of the elec- 
tron interaction up to second order in als- 

In Fig. 2 we have illustrated the form of the electron correlation by 
plotting the absolute value of C2D as a function of x in a log-log plot for 
a particular choice of the interaction strength and interaction length. The 
function is shown for different y-values, corresponding to points outside the 
edge, at the edge and inside the edge. For y < kpl B (which is 10 in the 
chosen units) the curves show a clear distinction between the short range 
bulk contribution, with a rapid fall-off with x, and the long range edge 
contribution with a slower, almost linear change with x. 

In Fig. 2, as well as in Figs. 3 and 4, the width W is measured in units of 
magnetic length, Vq gives the value of the dimensionless interaction strength 
Vq/A-kHvf and the inverse interaction length a is measured in units of inverse 
magnetic length. The correlation function C2D is measured in units of the 
bulk density l/27rZ^ and the coordinates x and y are measured in units of 
the magnetic length. 

5 Electron density and the asymptotic form of the 
correlation function 

The electron density is determined by the correlation function (58) for x = 
7] = 0. Since the function A(x) is proportional to x 2 for small x, we have 
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Figure 2: The absolute value of the electron correlation function C2d(x,U,v) as 
function of the relative coordinate x for different values of of the y coordinate. In all 
cases the y-coordinate is chosen to be the same for the two points of the correlation 
function (77 = 0). The value y = 10 corresponds to the edge value kpl% with the 
chosen units. 



A(0) = B(0) = 0, and the electron density is 



P2B(y) = ^MKl - a%C) y -±^] - erf[(l - a 2 /|C)^%61) 

with C = C(0). This constant is determined by the second derivative of 
A(x), and within the approximation applied here the explicit expression is 




V e~ z + AiThv F 



(V e- Z + 47rhv F ) 2 - e- 2a2w2 V$e- 2z 



(62) 



Thus, the electron interaction gives rise to a small change in the density, 
which is of second order in als- In the non-interacting case the density 
profile has the well-known form with constant density, corresponding to a 
full lowest Landau level, for y values between the edges at y = ikplg. We 
see from the expression (61) that the effect of the interaction is to widen 
the edge by introducing a slightly larger length scale Ib — > 'b/((1 — a 2 l B C), 
while the bulk density is left unchanged. For small a this length increases 
with a, but for sufficiently large a it is exponentially damped due to the 
presence of the factor exp(— 2a 2 W 2 ) in the expression for C. 
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a) 



b) 



Figure 3: The electron density as a function of the transverse coordinate y. Plot 
a) shows the density profile across the 2D system, with an almost constant density 
between the two edges. Plot b) shows the effect of the interaction, in the form of the 
difference between the density of the interacting and the non-interacting system for 
values of y close to the edge (at y — 10). The three curves correspond to three dif- 
ferent values of the interaction range a, with fixed value of the interaction strength 
Vq. The curve with largest deviation from the non- interacting case corresponds to 
aW = 0.6. 

The electron density is shown in Fig. 3 for a fixed interaction strength 
and variable interaction lengths. Fig. 3a shows the standard picture of a 
quantum Hall fluid, with a constant electron density in the interior of the 
2D channel, corresponding to integer filling, and with an abrupt fall in the 
density to at the edges. The width of the edge is essentially equal to the 
magnetic length Ib- For the chosen value of the interaction strength the 
effect of the interaction is not detectable in Fg. 3a, but is demonstrated 
in Fig. 3b where a magnified picture of the density is shown for y-values 
close to the edge. The graph shows deviations from the density of the 
non-interacting system for three different values of the interaction length. 
The maximal effect corresponds to an interaction length that approximately 
matches the width of the channel. 

We next consider the asymptotic form of the correlation function given 
by (58). We make use of the leading terms of the asymptotic expansion for 
the error function [11] 




(63) 



and the asymptotic expressions 





B{x) 



(7-1) — 
ax 
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C(x) « -(7-1)^2 ( 64 ) 



with xq as a constant. This gives 



C 2 D(x;y,rj) « - — 



— A+z 2 —A+z 2 



V 2 



2tt^ttI b x 4l B 

, {(y + k F lV) 2 , \ , (y-Ml) 2 
x(exp(- v ^ g; - iife F x) - exp(— ^ — g7 +ik F x)) 

B B 

(65) 

Near one of the edges, e.g. y w /cf^, this gives the asymptotic expression 

|C 2D (x;y,r,)| «K'^exp(-fc 7 J^)exp(-^-) (66) 

with as a new constant. The expression shows that the asymptotic 
x-dependence of the 2D electron correlation function is determined by the 
same exponent 7 as the ^-dependence of the ID (chiral) correlation function, 
(35). It also shows that the variation of the function in the y-direction has 
a gaussian form, with the function peaking at the edge of the system. 

In Fig. 4 we illustrate the asymptotic form of the correlation function 
and the transition from the small x behaviour, with the three curves showing 
the absolute value of the correlation function as a function of x for different 
values of the interaction strength. The constant values of the slope that 
are reached for large x are consistent with the corresponding values of the 
asymptotic parameter 7 in the three cases. However, we note that in the 
interacting cases a fairly large value of x is needed before the asymptotic, 
linear behaviour is reproduced. 



6 Conclusions 

We have examined in this paper the electron correlation function for a 2D 
quantum Hall system, which has the form of a linear channel of electrons in a 
strong perpendicular magnetic field. A simple model, which was introduced 
in Part I [1], has been used, where the electrons are confined by a harmonic 
oscillator potential and are subject to a gaussian two-particle interaction. 
An explicit expression for the 2D correlation function has been derived by 
use of a mapping of this system to a ID Luttinger liquid, a mapping that 
was discussed in detail in [1]. 

We have made use of the low-energy form of the theory in ID, where 
the electron correlation function can be determined by use of the bosoniza- 
tion technique, and have mapped this function back to 2D. In this way the 
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Figure 4: The absolute value of the correlation function C2D as a function of the 
coordinate x for fixed y = kpls (equal to 10 in the chosen units) and with 77 = 0. 
The interaction range is fixed, with als = 0.02, and the three curves correspond to 
three values of the interaction strength Vo, where one of them is the non- interacting 
case. The dashed curves are straight lines with slopes determined by the values of 
the asymptotic parameter 7 for the three cases. The plotted curves fit well the 
asymptotic lines for large x, but one may note that there are deviations for x 
values that are smaller than about 50 magnetic lengths, which is comparable to the 
range of the interaction potential. 



expression for the electron correlation function of the quantum Hall system 
has been found, which is valid within the low-energy approximation. The 
expression found in this way shows a natural distinction between a bulk 
contribution and an edge contribution to the correlation function. 

The method used here (and in Part I) is based on determining how the 
interactions in 2D modify the Luttinger parameters in ID, and how this in 
turn gives rise to modifications in the description of the 2D system. We find 
that the interactions only influence the edges of the system, since only the 
edge contributions to the correlation function are modified by the interac- 
tions. Thus, even if the particle interaction is long range and acts between 
the two edges of the system, there is only a negligible change in in the elec- 
tron density and correlation function in the interior of the quantum Hall 
channel. For the electron density even at the edge the effect is very small, 
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corresponding to a small widening of the edge. For the electron correlation 
function the main effect of the interactions is to change its asymptotic be- 
haviour. Thus, the fall-off with the distance x between the the two points 
of the correlation function is found to be determined by the same parame- 
ter value 7 that determines the asymptotic behaviour of the ID correlation 
function. However, for the specific values of the interaction strength and 
range that we have studied, the asymptotic form of the correlation function 
is reached only after a fairly long distance x, a distance that is comparable 
to the range of the interaction potential. 
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